function dlambda = ZermeloAdj(t, lambda, tu, u, ts, s)
% Backward integration
u = interp1(tu,u,t); 
s = interp1(ts,s,t); 

dlambda = [0; -lambda(1)];
